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ABSTRACT 

We report on the design, first observing season, and analysis of data from a new prototype millimeter- 
wave interferometer, MINT. MINT consists of four 145 GHz SIS mixers operating in double-sideband 
mode in a compact heterogeneous configuration. The signal band is subdivided by a monolithic 
channelizer, after which the correlations between antennas are performed digitally. The typical receiver 
sensitivity in a 2 GHz band is 1.4mK-y/s. The primary beams are 0.45° and 0.30° FWHM, with fringe 
spacing as small as 0.1°. MINT observed the cosmic microwave background (CMB) from Cerro Toco, 
in the Chilean Altiplano. The site quality at 145 GHz is good, with median nighttime atmospheric 
temperature of 9K at zenith (exclusive of the CMB). Repeated observations of Mars, Jupiter, and 
a telescope-mounted calibration source establish the phase and magnitude stability of the system. 
MINT is the first interferometer dedicated to CMB studies to operate above 50 GHz. The same type 
of system can be used to probe the Sunyaev-Zel'dovich effect in galaxy clusters near the SZ null 
at 217 GHz. We give the essential features of MINT and present an analysis of sideband-separated, 
digitally sampled data recorded by the array. Based on 215 hours of data taken in late 2001, we set 
an upper limit on the CMB anisotropy in a band of width A€ = 700 around I = 1540 of 6T < 105 fiK 
(95% conf.). Increased sensitivity can be achieved with more integration time, greater bandwidth, 
and more elements. 

Subject headings: cosmic microwave background — instrumentation: interferometers 



1. INTRODUCTION 

Measurements of the CMB at small angular scal es (i > 
1000) complement those from WMAP (Be nnett et al.l 
120031) and previo us experiments in cluding TOCO 
llMiller et all 119991) ■ B OOMERanG llMauskonf et all 
I2nnnt IRuM et al.l I200,'^ MAXIMA JFee et al.l I2001D 
DASI palverson e t al.' '200Z), VSA ('Grainge et al.l 
12003(1 ■ and ARCHEOPS l,Benoit et al. 2003), allowing 
yet more precise constraints on the standard cosmologi- 
cal model ()SDergel et al.ll2?ffl^ . Additionally, such mea- 
surements will enable us to understand the transition 
from the linear to the non-linear growth regimes for the 
formation of cosmic structure. One promising measure- 
ment technique is interferometry because of the ability 
to control syst ematic errors. Existing C MB interfer - 
omet ers, CAT llRobson et aLllT9 93'l. CBI llPadin et al 



ter using SIS (superconductor-insulator-superconductor) 
mixers at / > 100 GHz has the potential to measure the 
CMB anisotropy and SZ effects where the point source 
contamination is at a minimum fToffolatti et al."199'" 



19991) . 



DASI jLeitch et al.' '2002), "lAC (|Dickcr et al 
VSAp^atson et al. 2003) are based o n HEMT 
amplifiers fPospicszalski 1992; Pospicszal ski et al.ll99^ 

IPospieszalski igg?,) and operate at 30 GHz where the 

contamination from point sources can be large; MMIC 
technology at 90 GHz will enable a straightforward 
switch to higher frequencies. A compact interferome- 

^ Present address: NIST Quantum Electrical Metrology Divi- 
sion, 325 Broadway Mailcode 817.03, Boulder, CO 80305-3328. 

^ Present address: Department of Physics, University of Califor- 
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^ Present address: School of Medicine, University of Pennsylva- 
nia, Philadelphia, PA 19104. 

Department of Physics and Astronomy, University of British 
Columbia, Vancouver, BC, Canada V6T 1Z4. 

^ Present address: Center for Astrophysics and Space Astron- 
omy 389 UCB, University of Colorado, Boulder, CO 80309. 
NRC Fellow. 
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and where the SZ effect (Sunvaev & Zcl'dovich 197l_ 
119721) has a characteristic spectral null. However, the 
technical complications are considerable, as might be sur- 
mised by the existence of only a few such interferometers 
(e.g., SMA, OVRO, BIMA, IRAM, NRO/NMA). 

We have built and observed with a four element, het- 
erogeneous, compact, fully-digital, SIS-based D-band (~ 
145 GHz) interferometer designed as a prototype for a 
larger, more ambitious instrument. Novel elements of 
the interferometer include the channelizer, a custom mi- 
crowave integrated circuit; digital correlators for CMB 
measurements; and mixed antenna sizes. We describe 
key features of the interferometer, the observations in 
Chile, and the data analysis. 

2. the millimeter interferometer 

The Millimeter INTerferometer (MINT) is a compact 
instrument, approximately two meters on a side (with- 
out the ground shield) and cubical. All four receivers are 
mounted to a single rigid aluminum plate with fixed rela- 
tive alignment. The plate sits atop a movable aluminum 
cage, as shown in Figure ^ The channelizer, digital cor- 
relators, data computer, and a calibration noise source 
are also mounted on the cage. By moving the entire sig- 
nal path in unison, we avoid the potential phase changes 
that result from cable flexure and other moving parts. 

The telescope cage moves on an altitude-azimuth 
mount. A linear actuator drives the telescope in elevation 
over the range 60° to 100°, controllable to the accuracy 
of the 17-bit encoder (~ 0.005°). The azimuth drive has 
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Fig. 1. — MINT interferometer. An elevation over azimuth mount positions the cage that holds four Cassegrain antennas. 
A ground shield mounted on top of the cage moves with the receivers. External to the cage, and not shown, are a thermal 
regulation system that servos the coolant for the RF electronics and correlators to ±1 K, a position control computer, the helium 
compressors, and a 50 kW Katolight diesel generator to provide electrical power. The entire experiment fits inside one 20- foot 
shipping crate. 



a range of 270°, with an inaccessible zone centered to 
the north. The azimuth is also instrumented with a 17- 
bit encoder and can be positioned to ~ 0.015° accuracy. 
A brake on the azimuth bearing allows us to turn off 
the electrically noisy azimuth drive during fixed-azimuth 
observing. A ground shield made of aluminum hex cell 
in the shape of an inverted frustum extends 1.2 meters 
above the receiver mounting plate. This shield prevents 
the receivers from illuminating the ground, either di- 
rectly or through any number of reflections. 

2.1. Optics and Receiver Arrangement 

MINT uses antennas of two different sizes on its four 
receivers. A heterogeneous array is unusual in radio in- 
terferometry. It introduces complications such as com- 
plex effective beams and an "optical" phase that are un- 
desirable in most astronomy applications. However, the 
faint CMB anisotropy is a weak signal at the limit of de- 
tectability. We maximize the sensitive area in the aper- 
ture plane of the four-element prototype array, within 
structural constraints, by combining antennas of two dif- 
ferent sizes. Mounting the four antennas in a fixed con- 
figuration on a single plate also improves sensitivity by 
ensuring that the baseline separations in the {u, v) plane 
stay fixed as the sky rotates. Of course, deep observa- 
tions of a single (u, v) visibility come at the expense of 
mapping. 

An on-axis Cassegrain mirror system is mounted above 
each receiver. The inner two receivers have 30 cm di- 
ameter primary mirrors and 6.5 cm diameter hyperbolic 
secondaries. The equivalent focal length is 53.7 cm, giv- 
ing a focal ratio of //1. 8. The outer two receivers have 
Cassegrain optics of the same shape but larger by a fac- 
tor of 1.5 in each linear dimension. Rays from infinity 
focus near the entrance to a cold cylindrical feed horn* in- 
side the dewar. Identical horns with 17° FWHM beams 

* Custom Microwave, Longmont, CO 



feed all four antennas. The feeds are characterized in 
Miller et al. ( 2002). The mirrors were machined from 
aluminum on a computer-controlled lathe. A rigid tri- 
pod of thin (1.6 mm) legs made of G-10 composite sup- 
ports each secondary mirror (see Figure . A conical 
aluminum radiation shield with an opening half angle of 
5° surrounds each primary mirror to a height of approx- 
imately one diameter. The shields are rolled outward at 
the top with a radius of ~ 5A to limit diffraction and 
reduce crosstalk between receivers. 

The beam patterns of the two optical systems have 
been both measured and computed. Beam maps were 
measured using a chopped source (a 147 GHz Gunn diode 
with a D-band corrugated feed) and a diode power de- 
tector. The system allowed measurements over a dy- 
namic range of 35 dB. B eams were also com puted us- 
ing the DADRA program l|Barnes et al.ll200l .9 From a 
spherical wave expansion of the field from a corrugated 
feed, DADRA computes the currents on the subrefiector. 
Then, from the subrefiector, it computes the currents on 
the primary. The net pattern is the sum of the fields from 
the feed, subrefiector, and primary. We found it neces- 
sary to modify DADRA to include reflections from the 
primary back onto the secondary, because the MINT sec- 
ondary mirror has a larger diameter than the hole in the 
primary. The small primary hole is desirable, as it puts 
power onto the sky rather than backward into the de- 
war, but it must be included in the computation. Exten- 
sive measurements on the MINT and similar Cassegrain 
antennas (e.g.. Figure confirm the computed beam 
shapes. 

2.2. Receivers 

The receiver is shown schematically in FigureEl It uses 
the c lassic double -sideband superheterodyne configura- 
tion l|Kraus II 19861 e.g.). Other than the external optics, 

^ YRS Associates: Y. Rahmat-Samii, W. Imbriale, & V. 
Galindo-Israel 1995, rahmat@ee.ucla.edu 
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Fig. 2. — (a) Measured beam map of one of the small antennas. The dots are the measurement; the solid line is the DADRA 
computed beam, accounting for the defocus effect of measuring the beam at finite distance (36 meters). The data are normalized 
to the computed forward gain. Two-dimensional beam maps (not shown) confirm that the beams of both antennas are circularly 
symmetric. (6) Far-field beams, as computed by DADRA, for the small (solid) and large (dashed) antennas. 
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Fig. 3. — Schematic of the MINT RF chain, calibrator, and phase locking system. The noise calibrator is the top section; the 
sky signal flows from left to right in the center section; and the 145.1 GHz oscillator is produced and phase-locked in the bottom 
section. The dashed line surrounds cryogenically cooled components. 



the four receivers are essentially the same. Celestial ra- 
diation is focused by the Cassegrain optics onto a single 
corrugated feed. The feed couples to a polarizer that ac- 
cepts only right circular polarization. The output of the 
polarizer supports linear TEn mode in a circul ar waveg- 
uide (diameter = 0.173 cm). A transformer ijPadman I 
Il977|) converts it to TEio in WR-5.8. The ImW local 
oscillator^° (LO) at 145.1 GHz is coupled to the signal 
wave guide at —20 dB through a NRAO branch-line cou- 
pler ijKerr et al.1ll993,) . On the SIS junction, the LO 
(300 nW) and incident field (w 0.2 nW) are multiphed 
together. 

For computing the optical transmission into the feed 
horn, the ambient-temperature, 0.056 cm polypropylene 
vacuum window; the 25 K, 0.052 cm thick z-cut quartz 

Zax Millimeter Wave Corporation, San Dimas, CA 



IR blocker; and the 4K backshort are treated as a single 
multi-layer system. The spacing between the elements 
and the thickness of the window and quartz were selected 
to maximize the transmission near 145 GHz. 

Each SIS is cooled with a CTI-1020 CP refrigera- 
tor^^ modified with a third s tage f ollowing the design of 
iPlambeck. Thatte fc Svkes I ljl993|) . The modifications 
were made by Cryostar.-'^^ The resulting refrigerators 
are efficient; each has a cooling power of 0.05 W at 4K 
when driven by a 2.2 kW CTI 8200 compressor. The base 
temperature is approximately 3.7 K, but the temperature 
of the unloaded base varies by 400 mK peak-to-peak at 
25 Hz, the operating frequency of the fridge. Reducing 
the variation with a high heat capacity, high conductivity 

" CTI, Helix Technology Corp., Mansfield, MA 
Cryostar Associates, Yuma, AZ 
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combination of stainless steel and copper and judicious 
choice of thermal breaks, we were able to servo the cold 
head temperature to ±2 mK. This level of thermal stabil- 
ity was essential. Even with such temperature control, 
thermally driven dimensional changes in the cold head 
led to a small modulation of the reflection of the LO 
power from the vacuum window. The modulation could 
be monitored through ~ 0.5% variations in the SIS bias 
current. 

The SIS was tuned with a backshort coupled through 
vacuum-tight gears to a computer-controlled motor out- 
side of the dewar. The SIS is biased with a cryogenic 
NRAO bias-T, which outputs the RF signal. Follow- 
ing a 3.5-6.5 GHz cryogenic isolator, ^'^ a 4.1cm piece 
of stainless steel waveguide leads to an NRAO 3.5-6.5 
GHz, 5-8 K noise temperature, 33 dB C-Band HEMT- 
based amplifier at 14 K. The output of the C-band amp 
feeds through stainless steel coaxial cable to the room- 
temperature electronics where it is amplified by another 
52 dB. Half the signal is coupled to a total power monitor 
and the other half goes to the channelizer. With a lOK 
thermal load as input (similar to the load observed in the 
field) , the approximate level of the total power output is 
ImW. 

Each receiver is tuned while coupled to a cryogenic, 
thermally controllable blackbody load that matches the 
feed but is electrically and thermally isolated from it. A 
Lab VIEW computer program changes the load temper- 
ature, SIS bias voltage, and backshort position to find 
the optimal bias. In addition to the stable SIS temper- 
ature, the ambient electronics were stabilized to ±0.1 K. 
Repeated tests of the system over multiple cool-downs 
showed that the optimal tuning position was stable. The 
same tuning procedure worked without changes for all 
four receivers (A, B, C, D). Table^gives the noise tem- 
perature of the SIS system in the four bands. After initial 
tuning, the SIS bias was not retuned during the observ- 
ing season. 

To provide a continuous relative calibration, we built a 
138-142 GHz source that injects a broadband noise signal 
toward the subreflector of each telescope (Figure . A 
motor-controlled rotary mechanical switch ensured the 
"off" state of the source. 

The LOs for all four receivers are phase locked to 
a 100 MHz master clock, which also drives the digitiz- 
ers and correlators, the LOs in the channelizer and a 
12.1 GHz reference signal for the LO locking. In the 
ambient-temperature electronics associated with each 
receiver, the 145.1GHz LO is phased locked to the 
12.1 MHz tone using a circuit based on one given to us 
by the SMA project (Hunter ct al. 2002) at the Har- 
vard Smithsonian Center for Astrophysics. Each circuit 
is based on a type II phase-locked loop, servoing on the 
twelfth subharmonic of the 12.1 GHz-145.1 GHz mixing 
product. The LO frequency is controlled by a varactor, 
an approach that requires far less current than the more 
usual Gunn bias does. 

To summarize, the outputs of each receiver are a phase- 
locked 4-6 GHz IF signal and a total power monitor. The 
receivers were robust. They operated for a total of 30 

Passive Microwave Technology, Inc., Camarillo, CA and P & 
H Laboratories, Simi Valley, CA 

14 Pacific Millimeter Products, Golden, CO 



receiver- months without failure. The r eceivers are de- 
scribed in more detail bv iDoriesel 1)20021) . 

2.3. The Channelizers 

The channelizers subdivide the 4-6 GHz IF (interme- 
diate frequency) signal from each of the receivers into 
four 500 MHz bands, then downconvert each band to 
baseband. The channelizers are custom integrated mi- 
crowave circuits designed using Agilent EEsof Advanced 
Design System software. "'^^ Each circuit consists of a 
specialized layout of etched copper on 0.076 cm Duroid 
6200 substrate and drop- in components. Its only con- 
nectors (apart from power supplies) are on the input 
and output. The channelizers are housed in a custom 
aluminum enclosure that provides isolation between ad- 
jacent channels and from the outside environment. The 
cavities between the traces and enclosure are coated with 
adhesive-backed, rubberized microwave absorber to sup- 
press cavity modes. The channelizers for all receivers are 
mounted as a single, compact (6 cm by 30 cm by 30 cm), 
temperature-controlled unit. 

The splitting is accomplished by a four-way power 
splitter composed of planar cascaded Wilkinson power 
dividers followed by planar bandpass filters. The Wilkin- 
son splitters use a wid e-band, three-segmented design 
()Li. Li &: Bosisio1ll984|) . Each bandpass filter is a six- 
pole Che byshev type made from seven coup led microstrip 
sections (fMat thaei. Young fc JoneslllQSOD . The realized 
bandwidth varies from band to band, between 400 and 
600 MHz. There is no significant variation among the 
nominally identical bands in different channelizers. 

Each filter is followed by a commercial surface-mount 
mixer, which downconverts the signal to baseband (0- 
500 MHz), and a surface-mount amplifier. Two local os- 
cillator signals — one at 4.5 GHz and one at 5.5 GHz — are 
sufficient to mix the four separate bands to baseband. 
Thus half the bands are from the upper side-bands of 
the LOs and half are from the lower. The overall gain of 
the channelizer is ~ 12 dB across all four channels, with 
output power in the range +1 to +7dBm. We label the 
bands out of the channelizer according to their IF fre- 
quency: red, yellow, green, and blue. The red band is 
4.0 to 4.5 GHz in the IF, and thus nominally contains ce- 
lestial information from both sideban ds of th e SIS: 149.1 
to 149.6 GHz and 140.6 to 141.1 GHz. FTVMn l|2002f) gives 
further information about the channelizers. 

2.4. The Digital Correlators 

Each MINT digital correlator processes four 500 MHz- 
wide analog input signals, one from each receiver; all the 
red bands go to one correlator, for example. Incoming 
signals are digitized at 1 sample per nanosecond with 
2-bit resolution. The correlator produces a 16-lag cross- 
correlation function, which is read into a data computer 
twice per second. An offline Fourier transform of the cor- 
relation function produces a cross-power spectrum with 
62.5 MHz resolution. A total of four correlators are re- 
quired to process all sub-bands. The basic des ign of the 
corre lator was given to us by David Hawkins l|Hawkins I 

l2?iol . 

1^ Agilent Technologies, Palo Alto, CA 
1^ Mini-circuits, Brooklyn, NY 
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Table 1. Noise Temperature of Four Receivers Using NRAO SIS 

MIXERS''. 



Band Name IF Band (GHz) Ta (K) Tb (K) Tq (K) Td (K) 



Blue 


5.5- 


-6.0 


40 


33 


68 


29 


Green 


5.0- 


-5.5 


28 


22 


55 


31 


Yellow 


4.5- 


-5.0 


24 


19 


33 


27 


Red 


4.0- 


-4.5 


26 


19 


32 


24 



"Measured in the lab through the channclizcr with a regulated blackbody load. The IF band is nominal band of the SIS output; the 
bandwidth for each is close to 500 MHz. 



The electronics that compute the 16-lag cross- 
correlation functions between four analog signals are de- 
signed onto a single 38 cm by 38 cm, 8-layer, surface- 
mount circuit board. The four incoming signals are first 
processed by identical RF sections that each contain a 
variable attenuator, an amplifier, a low pass filter, and 
a voltage bias. The variable attenuators (M/A-COM 
model AT20-0107) have smallest steps of 0.5 dB. They 
are adjusted by data acquisition software to maintain the 
rms signal level into the digitizer near 0.28 V. This level 
optimizes signal-to-noise given our corre lator algorithm 
and the digitizer quantization at ±0.25 V ("Cooper ll97(Tl . 
Resetting the attenuators hourly was found in field tests 
to be sufficient. The low-pass filters (Mini-circuits SCLF- 
550) have a — 3dB point at 605 MHz and block unwanted 
noise (particularly leakage of the channclizcr LOs), which 
otherwise would be aliased into the signal. The voltage 
bias of —0.5 V is required by the digitizers. 

The RF pre-processor is followed by a commercial 
IGSPS, 6-bit digitizer (SPT7610).i^ The two most sig- 
nificant bits are passed to an ECL demultiplexer, which 
slows the signal rate to 62.5 MHz on a 32- bit bus. After 
conversion to LVCMOS, the 32-bit data streams from 
each receiver are fed into a single FPGA (Field Pro- 
grammable Gate Array, XCVIOOO),^^ which computes 
the correlation functions and communicates with a dig- 
ital I/O card in the data PC. The FPGA also keeps a 
running count of the digitized values for each receiver — 
the "digitizer histograms." These histograms, which are 
read out twice per second, proved very useful as moni- 
tors of the RF power reaching the digital correlator (see 

The algorithm for the correlator is built around a 
multiplier-adder kernel. The multiplier is implemented 
as a 4-bit lookup table, whose values are given by the 
4-level, "deleted middle products" correlator algorithm. 
This table has a quantization efficiency of e„ = 0.88 rel- 
ative to an ideal analog correlator (CooDcr 1970). 

The rack-mounted correlators for all receivers fit into a 
box 41 cm by 41 cm by 20 cm. Together the four consume 
800 W and are cooled with a thermally regulated fluid. 
The correlators worked as designed throughout the 2001 
campaign. More details are given bv .Tran (2002.'l . 

3. OBSERVATIONS 

3.1. Cerro Toco Site 

^"^ Signal Processing Technologies, Colorado Springs, CO 
^® Xilinx Inc., San Jose, CA 



Observations were made from an altitude of 5200 m 
on Cerro Toco^^ (longitude 67°47'09" West, latitude 
22°57'29" South), a mountain overlooking the high 
Chilean Altiplano near the ALMA and CBI sites. This 
was t he same site used for t he MAT/T OCO experiment 
(Mill er et all 11999: Torbet e t al.lll999fl . The mountain 
is accessible from San Pedro de Atacama, Chile, by four- 
wheel drive vehicle and can be reached from Princeton, 
New Jersey in 24 hours. 

The atmosphere at the Cerro Toco site was often excel- 
lent. During the 38 calendar days of the campaign (2001 
November 22 through December 29), twenty percent of 
potential observation nights were lost to overcast weather 
or snow. On the remaining usable nights, the effective 
atmospheric temperature at D-band was below 9K half 
the available time. The first and third quartiles were 7K 
and 13 K and the minimum observed temperature was 
5 K. The typical opacity in D-band was therefore only a 
few percent. These temperatures exclude the contribu- 
tion of the CMB and cover only the period from 23:30 to 
10:00 UT each night (local time is UT-3). Data totahng 
215 hours are used in the analysis presented here. 

Two radio communication systems and automated ob- 
serving software generally permitted remote operation 
of MINT from an apartmen t in the nea rby town of San 
Pedro, 35 km away (Miller eJaDl200l . One -way data 
transmission from the site (33 kbps data rate) allowed ob- 
servers to monitor the data in real time. A separate, two- 
way, low-bandwidth system employed packet modems to 
relay commands to the telescope control computer. 

3.2. Absolute Calibration 

The digital correlator measures the quantity p(t), the 
unitless correlation coefficient of two voltage signals as a 
function of lag r: 

{Va{t)V,{t + T)) 

where brackets represent averages over time t. The cor- 
relator out£uMs_conv2rted_^on^ following 
§8.3 of iThompson. Moran fc Swenson I (|1998j) . The next 
conversion, from p to correlated power or temperatures, 
requires knowledge of the overall system temperature of 
each receiver. 

The most robust power calibration comes from the to- 
tal power diodes, which measure power at the IF in a 
3 GHz band (wider than but inclusive of the signal band) . 

The Cerro Toco site of the Universidad Catolica de Chile was 
made available through the generosity of Prof. Hernan Quintana, 
Dept. of Astronomy and Astrophysics. 
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Small nonlinearities in the diodes are characterized in the 
lab. The linearized voltage is calibrated in the field to 
the input temperature, using a combination of hot/cold 
load tests and sky dips. The hot load is a large piece 
of Eccosorb at ambient temperature (T « 270 K); the 
cold load is the sky at zenith. By measuring the sky 
power at a range of elevations, we separate the sky tem- 
perature, the receiver temperature, and the receiver re- 
sponsivity (volts per Kelvin). The receiver temperatures 
and responsivities found in this way were consistent in 
three separate measurements taken at Cerro Toco over a 
three-week period. Assuming that responsivities remain 
constant, each diode serves as an independent monitor of 
the sky temperature. In practice, all four diodes found 
the same sky temperature to within 0.8 K rms, confirm- 
ing that no other important contribution to the measured 
power varied significantly during the campaign. 

The sky (atmosphere plus CMB) serves as a variable- 
temperature load, given the Eccosorb and sky-dip cali- 
brations. The digitizer histograms give the power in each 
0.5GHz-wide sub-band. Together, the power and sky 
temperature permit an estimation of the receiver temper- 
ature and responsivity of each sub-band. Although the 
responsivity is found to change in one or a few discrete 
jumps during the campaign for some bands, the noise 
temperatures do not change significantly. Thus the total 
system temperature in each channel can be determined 
at any time simply by adding this constant receiver noise 
temperature to the varying sky temperature. The jumps, 
of approximately 10-20% in magnitude, occur only be- 
tween nights and are attributed to problems in signal 
cables or their connectors. 

This calibration method assumes that several key re- 
ceiver properties stay constant: the responsivity of the 
total power diodes; the receiver temperature at the 
diodes; and the receiver temperature in each 0.5 GHz- 
wide sub-band. The redundant measurer nents of the sk y 
temperature confirm these assumptions (jPoriese 

Observations of Mars and Jupiter are used to estab- 
lish the forward gain (including optical efficiency) of the 
receivers. The gain is required to convert between bright- 
ness temperature and the fluctuations in the CMB phys- 
ical temperature. Planetary observations are described 
in Section^ 

3.3. Relative Calibration and Phase Stability 

The phase and amplitude stability of the interferome- 
ter are monitored using a switchable noise source. The 
noise is produced by a -t-20dB ENR Q-band source^" fol- 
lowed by a series of amplifiers providing a gain of -f 70 dB. 
This signal enters an isolator, a frequency tripler, and 
another isolator (Figure A D-band rotary waveguide 
switch (controlled by a stepper motor) allows the noise 
to be turned on or off. The signal is then split four ways 
by a 180° magic-T hybrid coupler and a pair of 90° short 
slot hybrid couplers. The resulting four signals are highly 
correlated, except that each is phase-delayed by 90° rel- 
ative to the next. All components up to and including 
the power splitters are housed in a temperature-regulated 
box. Each signal travels through one meter of coin-silver 
D-band waveguide, where it is coupled into the telescope 
optics by an open-ended waveguide illuminating a tiny 

Noise Com, Parsippany, NJ 



"pickoff" mirror adjacent to the dewar window. Some 
power reflects up to the secondary and back into the de- 
war. The coupling efficiency is not measured but is found 
to be stable. The noise source is sharply cut off above 
145 GHz by the bandpass of the Q-band amplifier, with 
the consequence that the noise contains D-band power 
only in the MINT lower sideband (139 - 141 GHz). This 
property is helpful in verifying the sideband separation 
procedure f H6.4|l . 

During normal CMB observation, the noise source is 
switched on for 6-8 seconds out of every 450. A longer 
calibration of 30 seconds is performed every hour. The 
noise source produces power equivalent to 5-10 K in each 
receiver, so the correlated signal is observed with high 
signal-to-noise in a single 2-second integration period. 

The magnitude of the noise signal recorded by the in- 
terferometer was stable throughout the campaign. On 
most nights, the fractional range of amplitudes was 3% 
(rms) or less in all bands and all baselines. Over the 
entire campaign, the rms range was 6%. Strong correla- 
tions in the amplitudes among all bands and all receivers 
suggest that the main source of amplitude variation was 
the noise source itself. For example, some variation is due 
to the waveguide switch failing to find the same resting 
spot on each repetition. 

The phase stability of the instrument — a critical re- 
quirement for any interferometer — was also measured us- 
ing the noise source and was found to be sufhcient. Con- 
sidered over the campaign, the measured phase variation 
over time and across frequency bins would cause average 
signal loss of 6% in most baselines (corresponding to an 
rms phase variation of 20°). The D-C and A-B baselines 
performed better than the others, with average loss of 
less than 3% and (prms = 12°. We attribute this differ- 
ence to the power-splitting scheme in the noise source 
(A and B signals come from the same 90° hybrid) , which 
permits a relative phase between the A/B signals and 
the C/D ones. This observation suggests that the cali- 
bration source itself is responsible for the main part of 
the observed phase variation over time. Therefore, the 
interferometer is more phase-stable than could be mea- 
sured with this source (a comment that also applies to 
amplitude stability), and the 12° rms phase variation is 
an upper limit on the stability of the interferometer itself. 

The instrument and sky noise is observed to integrate 
down with rms signal values scaling as t~^/^ for averaging 
times t of up to 2000 seconds in a single night. Longer 
integrations were too few to make significant comparisons 
with the expected scaling law. 

4. OBSERVATIONS OF THE SUN, MARS AND JUPITER 

Repeated observations of the Sun were used to deter- 
mine the telescope alignment throughout the Chile cam- 
paign. The receivers were temporarily covered with stan- 
dard corrugated cardboard 2 mm thick to block the Sun's 
optical and infrared radiation. The pointing was derived 
from separate maps of each receiver's total power, not 
from interferometry. Common-mode pointing errors can 
be corrected off line and were determined to approxi- 
mately 0.025° accuracy in elevation and in azimuth. Fol- 
lowing an initial alignment of the receivers to each other 
(which also involved solar observations), they were found 
to be relatively aligned to within ~ 0.02°. This angle cor- 
responds to approximately 1/6 of the fringe spacing in 
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the longest baseline. We noted no significant shifts dur- 
ing the campaign. At this level of accuracy, the flatness 
and placement of the cardboard shields appear to limit 
the measurements. 

During the Chile campaign. Mars was visible nightly 
in the late afternoon and early evening. When weather 
and maintenance permitted, MINT observed Mars by re- 
peatedly pointing the telescope ahead of Mars' position 
and allowing the planet to drift through the beam for 1-2 
minutes. Approximately twenty hours of data were taken 
on or near Mars. Averaged over the season, the planet's 
angul ar diameter was 7". Using the 171 GHz temper- 
ature l|Goldin et al.lll997fl of 198 K, we expect Mars to 
produce antenna temperature in the small-small, mixed- 
size, and large-large baselines of only 1.9 mK, 2.9 mK, 
and 4.3 mK (the three baselines have forward gains of 
-K51.3dB, -h53.1dB, and -1-54.80 dB, according to the 
lossless optical model, and solid angles of 93 /xsr, 62 /zsr, 
and 42 /isr). Jupiter was observed only at the end of the 
campaign. 

The Mars data are used to measure the phase versus 
frequency profile of the instrument. Since Mars appears 
as an unresolved point source to MINT, its visibility 
phase should be a simple function of its position rela- 
tive to the phase center of a given baseline. In particular, 
the phase should be nearly independent of the base-band 
frequency (effects of wavelength variation are small due 
to the small fractional bandwidth at D band). With 
the phase corrected for the position of Mars, a fiducial 
phase (f>{i') is determined by averaging the complex visi- 
bility over many Mars observations. All data, including 
planetary calibrations and CMB data, are adjusted by 
subtracting the fiducial phase. After making this adjust- 
ment, it is possible to average complex visibility across 
all frequencies. 

The visibilities measured on Mars match the impor- 
tant features of a point source. The visibility magnitude 
is approximately Gaussian (as a function of pointing an- 
gle away from the planet) with a width as expected from 
the primary beam widths. The visibility phase is also 
consistent with the expected 2ttx • u fringes. The phase 
of the setting Mars is stable over the campaign at the 
level of 5°-10° rms, depending on the baseline. Even 
in the worst baseline, the computed signal loss due to 
the phase variation is only 3%. The phase is somewhat 
less stable during observations of Mars rising (10°-20° 
rms, or 2%-10% loss), but these data are taken during 
daylight hours and are therefore less indicative of the 
instrument stability during its nighttime CMB observa- 
tions. 

Mars' visibility magnitude, though stable during the 
campaign, does not match the expected magnitude. For 
all baselines involving receiver A, the magnitude is only 
33% of the expected value; for other baselines, the obser- 
vation is 75% of the expected value. Extensive tests to 
explain the deficit were performed after MINT returned 
to Princeton in 2002. The tests involved direct observa- 
tion of Jupiter, a high signal-to-noise source, and con- 
firmed the signal loss. By switching equipment between 
the receivers, channelizers, and correlators, the deficit 
was localized to the poor performance of the circular-to- 
linear polarizers at the base of the feeds. The degradation 
is clearly evident only while observing a known celestial 
source interferomctrically (which was not possible prior 



to deployment). This loss of coherent signal is corrected 
at the end of the analysis pipeline by scaling up all aver- 
age and rms values. 

Jupiter was too low in the northern sky during late 
2001 to be visible to MINT in its normal configuration. 
After the CMB campaign ended, the ground shield was 
removed and replaced by a flat mirror that covered all 
receivers. This mirror redirected the beams of each re- 
ceiver by 76° in altitude, allowing observation of Jupiter. 
The planet's larger angular diameter (47") makes the 
expected signal much larger than the Mars signal and 
easily discerned in a single 2-second data frame. We 
assume a brightness temperature of 172 K for Jupiter, 
given by WMAP measurements at 90 GHz l|Page et al.l 
'2003). The Jupiter observations (taken on consecutive 
mornings 2002 January 2-4) confirm all the main results 
of the lower signal-to-noise Mars data, including phase 
and amplitude stability, beam sizes, and the reduction 
in signal (particularly in receiver A) caused by the po- 
larizers. 

5. OBSERVING THE ANISOTROPY 

5.1. Sky Scan Strategy 

An experiment's sensitivity to CMB anisotropy de- 
pends on how the observing time is divided among fields 
on the sky. The statistically optimal division requires 
observing each field until a signal-to-noise ratio of one 
is achieved, then moving on to the next. For MINT, 
this meant maximizing the observing time on each field. 
We used a quasi-tracking strategy which combines the 
stability advantages of the sky drifting over a stationary 
telescope with increased integration time beyond a nor- 
mal drift-scan. The telescope was repointed eight times 
per hour to observe the same strip of sky (1/8 hour wide 
in right ascension) repeatedly. Pointings were made at 
a fixed azimuth angle and a range of near-zenith eleva- 
tions. Although fixed points in the sky curve across a 
range of azimuth, this effect was small for the limited 
range of elevations used. 

Table 12 gives the location of the center of each strip 
used in the current analysis, both in ecliptic and in galac- 
tic coordinates. Although the fields near RA=7^ and 8^ 
are within 10° of the galactic plane, they are located 
far from the galactic center and outside other major 
components. Additionally, the spectrum of diffuse back- 
grounds such as interstellar dust falls approximately as 
(Gauticr et al. 1992). The galactic contribution in 
the MINT fields at multipoles i > 1000 is negfigible. 

The effect of scanning continuous strips is to acquire a 
series of complex visibilities with signal highly correlated 
between neighbors. Simulations show that sky strips w = 
1.875° long observed with a Gaussian primary beam of 
width (Tb = 0.19° have the same statistical power for 
measuring sky variance as rispot = 4.5 separate spots. 
The general relation for all three MINT beam sizes is 
found to be 



5.2. Data Selection 

The data used for CMB analysis are selected based on 
time of day, sky temperature, and nominal receiver op- 
eration. The noise source calibrations f H3.3|) show that 
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Table 2. Median right ascension'' and galactic coordinates for CMB 
fields with at least 9 hours observing time*'. 



RA (h,m 


.s) 


6 (deg) 


I (dog) 


1 3 


42 


—84.47 


161.59 


2 3 


5 


-72.68 


202.49 


3 2 


40 


-59.54 


212.27 


4 2 


34 


-46.34 


218.46 


5 2 


26 


-33.12 


223.82 


6 2 


16 


-20.24 


229.17 


7 2 


8 


-7.68 


234.97 


8 2 


6 


4.37 


241.77 


9 2 


1 


15.52 


249.94 


10 1 


47 


25.41 


260.21 


11 1 


39 


33.31 


273.18 


12 1 


38 


38.34 


289.07 



''The declination of each strip is 23° 2' 14" south. 

''The nominal fields are spaced evenly in right ascension. The coordinates given here are medians of the data surviving all cuts. 



the system phase response is stable from 23:35 UT un- 
til 10:05 UT. These times correspond approximately to 
half an hour after sunset and twenty minutes after sun- 
rise. Outside these times, the phase variations are likely 
due to thermal changes in the noise calibration waveg- 
uide (which is insulated but not thermally regulated). 
Although we believe that the instrument is more stable 
than the calibrator, there is no way to be certain, and 
data are rejected outside this window of 10.5 hours per 
night. 

In addition, the data are cut when the sky brightness 
in D-band (including CMB) exceeds an effective temper- 
ature of 35 K. At higher temperatures, the linear re- 
lationship between sky temperature and measured rms 
visibilities breaks down. The sky brightness cut removes 
3% of the otherwise valid data, which due to their high 
noise temperature would have only 1% of the statistical 
weight. 

Instrument criteria include tests for locked PLLs, func- 
tional 4K refrigeration, and properly synchronized cor- 
relator readout. The total data set after all cuts is 215 
hours long, representing 54% of the available nighttime 
hours. The weather and the instrument are responsible 
for approximately equal shares of the lost time. 



6. ANALYSIS 



6.1. 



Expected CMB Signal for a Heterogeneous 
Interferometer 

MINT has both 30 cm- and 45 cm-diameter antennas, a 
configuration which produces high collecting area for the 
medium and longest baselines, where the CMB signal is 
small, while permitting the shortest baseline to be only 
160 A. Put another way, the heterogeneous array fills a 
larger fraction of the antenna plane than would an array 
of small antennas alone. 

Given a CMB power spectrum, the expected vari- 
ance of the interferomet er visibility is derived by 
iHobson. Lasenbv. fc Jonesi l)199,5D and iWhite et al.l 
l)1999D for the case of identical receivers. We consider 
here the calculation for a heterogeneous array, which dif- 
fers in that the effective beam and its Fourier transform 
are complex. 

We work with a modified "temperature visibility" 
quantity, which corresponds to a brightness tempera- 



ture averaged over the beam solid angle, rather than the 
usual flux density. This definition amounts to a re-scaling 
of the usual visibility by a constant factor / {2knB), 
where fis is the beam solid angle. We make this non- 
standard choice even though MINT's power levels are 
ultimately calibrated to the flux density received from 
point sources (Mars and Jupiter) . We have found it con- 
venient to work in sky brightness temperature, a quantity 
appropriate to a beam-filling source. In particular, this 
choice permits quick and transparent comparison with 
non-interferometric CMB experiments and their sensi- 
tivities. We use the label Vr to designate the modified 
visibility definition; all formulas for Vr can be multiplied 
by 2kil. b/X^ to yield the conventional visibility. 

Ignoring the frequency dependence of the CMB bright- 
ness across MINT's narrow band and making the flat-sky 
approximation for small fields of view, the complex tem- 
perature visibility Vr at location u in the u — v plane 
is 

Vt{u) = / —A{x)ATB{x)e+^~, (1) 

where x is the two-dimensional position in the plane of 
the sky (conjugate to u), A is the primary power beam 
of the antennas (normalized to one at its peak), ^Ib = 
J dx A{x) is the beam solid angle, and ATb is the CMB 
brightness temperature fluctuation (ATb — rjj^jAT). 
The Rayleigh- Jeans factor rjnj converts small physical 
temperature fluctuations AT to brightness temperatures: 

VRj = ~ 2fc 9r ~ (e^ - 1)2 ' 

where x = hv/(kT). At 145 GHz and for T = 2.73K, 
VFtj — 0.59. Note that the definition of temperature 
visibility (equation ^| leads to a magnitude of \Vt\ — 
Tb ri/f2_B for an unresolved source of brightness temper- 
ature Tb and solid angle Vl <C ^b- This fact connects 
the CMB observations to the absolute calibration off of 
Mars. 

Apart from constant factors, the visibility is the 
Fourier transform of the product of A and ATb ■ By the 
convolution theorem, it is also the convolution of their 
transforms: 

Vt{u)^'^{A®^){u). 
I Lb 
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Using this expression for temperature visibihty, its ex- 
pected covariance matrix has diagonal elements of 



--.{V^{u)Vt{u)) 

( IBlY 



(2) 



dv\A{u^v)\'T^,„,S{v) (3) 



using the unitless power spectrum, S{v), defined by 

Note that the Fourier power spectrum 5* depends only on 
the magnitude of v and approximately equals the spher- 
ical harmonic spectrum Ci for large i (if the spectrum is 
smooth enough): S{v) « Ci^2ttv 

Non-diagonal elements of the theory covariance ma- 
trix are of two types. Two observations made by a single 
baseline but separated by angle A on the sky are accom- 
modated with a factor exp(27riA • v) in the integrand of 
Equation|21 Observations made by distinct baselines are 
also correlated. In the extreme case of parallel baselines 
(such as DA-BC and AC-DB in MINT), the correlation 
is complete. The MINT baselines of equal length but 
different orientations (such as DB and DC) are weakly 
correlated. The ratio of the magnitude of the covariance 
for DB-DC data to that of the DB-DB data is 18% for 
data from a single field (A = 0), falling off quickly with 
increasing field separation A. The correlations between 
the short or the long baselines (DA or CB) and any other 
are negligibly small. 

Analysis of the data is more easily done using the real 
and imaginary visibilities separately, leading to the final 
covariance elements 



^RR 



and 



C 



RI 



2n% 



2nl 



dv \A{u - v)\'^S{v) cos(27rA • v) 

(4) 



sin(27rA • v). 



(5) 

The other combinations obey Cfj = Cj^j^ and Cj^ = 
— C^j. The new factor of 2 accounts for the fact that 
half of the sky variance appears in the real component of 
V and half in the imaginary. The two formulas are even 
and odd with respect to exchange of the two indexes 
(or equivalently, under reflection of A). This result can 
be generalized to give the covariance between distinct 
baselines. 

In words, the covariance for A = is the sky power 
spectrum S(v) averaged in two dimensions with a weight- 
ing function |A(m — t>)p. Recall that A{x) is a single re- 
ceiver's power beam, so A{u — v) is the autocorrelation 
of the voltage pattern in the antenna plane, offset from 
the origin by a distance u corresponding to the baseline 
separation in wavelengths. 

The complication of mixed antennas can be handled by 
extending the formulation of the antenna power reception 
pattern fThomDson, Moran & Swcnson 1998). Starting 
instead with the voltage pattern gj{x) of antenna j for 
radiation from angle x, the effective power beam Ae for 
baseline j-k is 



With this change, all formulas and statements above 
about A and A hold, except that A{u — v) becomes 
the cross-correlation of the two antenna's voltage pat- 
terns in the aperture plane. The resulting beam Ae is, 
in general, complex. The overall phase is unimportant; 
it merely adds to any other phase offsets in the sys- 
tem, which are measured by observing planetary point 
sources. Phase variation across the beam, however, adds 
an imaginary part to the beam A and its Fourier trans- 
form A. In mint's mixed-antenna baselines, the phase 
of the effective beam is nearly constant well past the half- 
power point of the main lobe, so the imaginary part of A 
is of much smaller magnitude than the real part. Even 
so, the imaginary part is included in \A\'^ when numeri- 
cally integrating equations01andEl increasing the theory 
covariance by a few percent. 

6.2. Approximate Analytic Expression for the Expected 

Signal 

The theory covariance integrals can be approximated 
by taking the effective beams to be Gaussian and the 
power spectrum to be scale-invariant. The result is use- 
ful in illustrating how the CMB signal depends on in- 
strument parameters. We assume a Gaussian beam with 
a beam solid angle fie = ^ira'^. The ratio of the beam 
solid angles for the idealized Gaussian to the actual beam 
is akin to a main beam efficiency; call it eg = 2710"^ /VLb- 
The transform A of the normalized beam is then 



Let 



A{u) = 27r(T^ e 



S{v) 



27r7i2 V 



5T^ 



so that 5T'^ is a flat-bandpower temperature variance. 
Performing the angular part of the integral over dv 
(Equation |3Jl yields 



_ iVRJ eg STY 

'-^RR — 7^ 



dv 



e-(2--)'(«'+-')/o(87rV2w«), 



where In is the modified Bessel function of t he first kind 
l|Dieguez llT999t iHobson fc MaisingeTI 120021) . This can 
be simplified with the asymptotic expansion Iq{x) ~ 
e^ I \/2ttx for large x. Defining U = 2TTau and V = 2TTav, 



^RR — 



iVRj eg STy 
2C/2 



dV 



1 



U 

An \V 



3/2 



,-{u~vy 



Mx) = gjix)gl{x). 



Note that rp — A/ (27rcr) is the characteristic radius of the 
Gaussian power reception pattern in the antenna plane. 
Therefore, U = Xu/rp > 2 for any real interferometer, as 
the separation Xu must be larger than the antenna diam- 
eter. The integrand above is small except near V ~ U, so 
the divergence as ^ can be ignored (it corresponds 
to an unphysical divergence in the CMB power spectrum 
at large angular scales, where the flat-sky approximation 
breaks down as well). For U ~ 3.5, the integral evalu- 
ates numerically to 0.54, approaching 0.500 in the limit 
of large U (large separation). Replacing the Gaussian 
beam width cr with the b eam full width at half-maximum 
power OpwHM = \/8 In 2 ct (which is more readily gener- 
alized to non-Gaussian beams), the simplified estimate 
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for the covariance matrix is 

C]ij,^ 0.27 el ST' (6) 

^ 0.038 tjIj el STy{u9FWHMf- (7) 

Equation \7\ can be used to estimate the magnitude of 
the CMB signal from any single interferometer baseline, 
given basic parameters such as baseline length and the 
beam width and efficiency. We expect that the "Gaussian 
efficiency" eg can generally be replaced with the main 
beam efficiency with little loss of accuracy. 

6.3. Expected Signal-to-Noise on the Power Spectrum 

A numerical evaluation of equations 0] and [3 gives a 
more accurate figure for the expected MINT variance per 
fxK^ of flat bandpower (Table OJ. The two-dimensional 
Fourier transform of the beam, A, is simplified by noting 
that the antennas and beams are symmetric about the 
optical axis. This symmetry reduces the problem to the 
one-dimensional Hankel transform of the radial function. 
The numerical result assumes the beam shapes computed 
by physical optics and confirmed by direct measurement 
(Figure Ell. 

The noise covariance follows straightforwardly 

from the system temperature and the Dicke equa- 
tion |bicke , T9461 . The covariance of the real or imagi- 
nary output of a complex correlating interferometer is 

^ (ATI) = 

where Tgyg is the system noise temperature, is the ef- 
ficiency factor of the digital quantization scheme, Av is 
the IF bandwidth, and t is the effective observing time 
(see for example Jhompson. Moran & Swcnson (1998) 
eq. 6.42). Note that "for MINT, r is only half the 
elapsed time, because the phase-switching scheme effec- 
tively time-shares the correlator between measuring the 
real and imaginary parts of the correlation. Using a 
typical Tgys = 45K (including atmosphere and CMB), 
eg = 0.88, Av — 1.9 GHz, and 5 hours of observing time 
(r = 9000 s) on a CMB field, the noise covariance is 

Using a set of iV ^ 1 Gaussian deviates to estimate 
the variance of the set gives a fractional uncertainty of 
y'2/N. The number of samples N is twice the number 
of fields observed, n/, because the complex correlation at 
each field gives two independent samples. The variance 
C'^ + of the data can be estimated to a fractional 
uncertainty of \/ ^JnJ, so 

The two terms are the sample and noise variances. The 
quantity of interest, i5T, is an rms rather than a variance; 
thus if the error is small, ST has half the fractional uncer- 
tainty of C'^ . Consider the case of MINT observations, 
covering Uf = 50 fields. If sample and noise variance 
contributed equally to the final error (which is statisti- 
cally optimal for a fixed total observing time), then the 
fractional error expected on the CMB signal ST would be 
l/\/50 ~ 15%. Achieving this level of uncertainty would 
require 76 hours of observing per field, approximately 12 
times the time actually available in the 2001 campaign. 



6.4. Sideband Separation 

The SIS mixers in MINT arc inherently double- 
sideband devices, mixing sky signals in both ranges 139- 
141 GHz and 149-151 GHz down to a single IF band 
at 4-6 GHz. By switching a 90° phase shift in and 
out of the 145.1GHz local oscillator tone every half- 
second, we separate the signals from the two sidebands. 
The separation is possible because an LO phase shift 
enters with opposite signs in the IF signals from the 
two sidebands. If C and C are the correlation func- 
tions for unshifted and shifted data, then (C -|- iC) and 
(C — iC) isolate the contributions from lower and upper 
sidebands. This feature is preserved through the chan- 
nelizer, where the signals are mixed do wn a second time. 
iThompson. Moran fc Swenson I l|1998f) describe the tech- 
nique in §6.1. We have shown that their result generalizes 
to the cases of unequal power in the two sidebands and 
of lower sideband conversion at the second mixing stage. 

The magnitude and phase of the cosmic signal should 
be nearly the same in both MINT sidebands. There- 
fore, it is not strictly necessary for CMB observation to 
distinguish the signals in the two bands. We could main- 
tain a constant LO phase, effectively measuring only one 
component of the complex sky visibility. Instead, MINT 
determines the full complex visibility with half the ob- 
serving time per component. The MINT strategy was 
chosen to characterize the instrument completely and to 
guard against possible drifts in the instrumental phase. 
Correcting such phase drifts — if measured by the calibra- 
tion source (§ I3.3|l — requires the complex visibility. 

The sideband separation technique is proved in prac- 
tice by observations of the calibration noise source. The 
noise source produces power only in the lower sideband of 
the MINT receivers. The measured cross-power spectra 
of the calibrator shows no power in the upper sideband 
except for small leakage near the band edges, a result of 
the limited number of time lags available in the correla- 
tion function. 

7. ANISOTROPY RESULTS 

7.1. Data Reduction 

Analysis of the interferometer data is performed in 
Fourier (temporal frequency) space. First, the correla- 
tion function p(t) is converted to brightness temperature 
(as described in ti3.2l and 2J ■ The conversion incorpo- 
rates the results of sky dips and the absolute calibra- 
tion determined from Mars observations. Using the side- 
band separation procedure, data with and without the 
90° phase shift in the LO are combined to form a com- 
plex correlation function. The discrete Fourier transform 
of the 16-lag correlation gives the cross-power spectrum 
of one pair of antennas. Interpretation of the cross-power 
is complicated by the fact that both the lowest and high- 
est frequency bins combine signal from upper and lower 
sidebands. This complication is not a problem, because 
MINT does not currently make use of the spectral infor- 
mation in analysis of the CMB data. 

The "fiducial phase" 4>{vj) found from many Mars 
observations (SQ is subtracted from each data point. 

The lower sideband power exceeds the upper by 7% for CMB 
fluctuations, owing to the frequency dependence in the conversion 
from physical to brightness temperatures. 
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Table 3. MINT Baselines, Bandpower, Conversion Factors, and Re- 
sults OF 2001 Campaign.'' 



Baseline 


u 


Beam solid 
angle Slg 
{^lsr) 


Beam Conversion factor, 
FWHM Gaussian 
dFWHM (deg) (10-3 mKVmK^) 


Conversion factor. 
Numerical 

(10-3 mkV/^k^) 


Effective £ for 
flat ST 


Upper limit on &T 
(pK, 95% C.L.) 


Short 


169 


93 


0.445 4.3 


4.3 


980 


232 


Mixed'' 


260 


62 


0.348 2.5 


2.4 


1540 


105 


Long 


492 


42 


0.297 1.1 


1.0 


3030 


323 



''The two conversion factors give the ratio of signal variance (in brightness temperature) to the CMB bandpower C^/5T^. For example, 
a bandpower of 1000 /xK^ would produce a signal of 4.3 /iK^ in the short baseline. Gaussian and Numerical columns refer respectively to 
the simplified estimate (Equation |2|in the text) and to the full numerical integration of Equation loused for the analysis. The conversions 
include several factors. The Rayleigh- Jeans factor squared is ~ 1/3. Signals separated into real and imaginary parts each contain only half 
of the anisotropy power. Use of a realistic beam instead of a "top hat" introduces a further factor of two in signal. Finally, the largest 
contribution, a factor of 15-60, comes from incomplete sampling of the u-v plane. The effective £ is the mean value probed if the power 
spectrum is flat. The last column gives the upper limits found by MINT on the CMB anisotropy. 

''The conversion factors are per baseline. MINT has four baselines of the intermediate length u = 260. 



The subtraction accounts for the constant phase varia- 
tion across frequency of components such as imperfectly 
matched waveguides and amplifiers. The visibility mag- 
nitude is also divided by a template (based on plane- 
tary observations) to correct for gain variation across the 
signal band. With 16 lags and four separate 500-MHz 
sub-bands, the template for each baseline has 64 com- 
plex elements, which are determined once for the whole 
campaign by averaging the Mars results. Phase- and 
magnitude-corrected data are now directly comparable, 
and are therefore averaged across all 64 frequencies. 

After all cuts, there are 390,000 two-second data 
records. Each record is a complex visibility (in brightness 
temperature) for each of the six baselines. The data show 
non-zero averages over the campaign of approximately 5- 
20 jiK. These offsets do not vary significantly with time 
or with the observing elevation, so the six complex off- 
sets are subtracted from the data records. Cross-talk 
between signal channels is the most likely cause of the 
small offsets. 

CMB data are combined into bins of width 0.1° in right 
ascension (the bins are smaller in actual sky angle by a 
declination factor of cos (5 = 0.921). All data are taken 
at fixed declination, so sky binning is required only in 
one dimension. The bin size is chosen to be the largest 
bin with nearly complete correlation of signal visibilities 
across its width. The visibility of a CMB signal varies 
across the bin with the usual 27rM • x phase, but this po- 
sition dependence is removed by "phase-centering" the 
visibility as if it had been observed at the center of the 
bin. Given a width of 0.1°, the signal lost due to av- 
eraging is only 1.5% in the longest baseline and 1.2% 
and 0.8% in the others, a loss which is accounted for by 
scaling the final results. The signal reduction would be 
closer to 10% for most baselines if the phase-centering 
correction were not made (the exception being the two 
baselines with u oriented north-south). 

Combining the data within a bin entails finding the 
weighted average and the uncertain t y on the mean, as 
given in IHobson. Lasenbv. fc Jones I l|1995() . Each data 
point has a noise variance given by the 15-minute average 
of the visibility variance, scaled by the instantaneous sys- 
tem temperature to account for variation in the sky tem- 
perature on faster time scales. The rms estimated from 
this procedure matches the rms of the phase-centered 



data to better than 2% for all baselines, when taken over 
the entire campaign. For the final likelihood analysis, 
the estimated "instantaneous rms" is scaled to eliminate 
the small mismatch. 

Data from the identical baselines DC and BA could 
be averaged together at this point, but the large loss in 
Receiver A makes the BA data subject to unknown sys- 
tematics and far less sensitive than the DC data; the 
AC data is likewise less sensitive than the BD data. We 
therefore discard the results from the BA and AC base- 
lines. 

7.2. Results 

The time-averaged visibilities at each bin are our 
best estimate of the complex sky visibility on a 
0.1° R.A. spacing. A maximu m likelihood analysis 
(e.g.,fe ond. Jaffe fc Knox I l|2000D : iHobson fc Maisinseri 
()2002j) ) is used to find the fiat bandpower in three dis- 
tinct bands, corresponding to the MINT baseline lengths. 
The signal covariance matrix, given by equations 21 and 
|31 is computed numerically for the appropriate power 
spectrum {S{v) cx w^^) for all six baselines, and for all 
possible separations A from zero to 3° in 0.1° steps. For 
this computation, the beam transform A{u) comes from 
the physical optics code. The DB-DC correlation of 18% 
is included in the covariance matrix, but the small co- 
variances between baselines of unequal length (< 1% in 
all cases) are ignored. Because the scan strategy concen- 
trated on 2° strips spaced every 15°, the signal covariance 
is block diagonal with approximately 18 separate blocks 
(of the 24 possible blocks, six correspond to areas of the 
sky high during the afternoon, which were not observed). 
The noise covariance matrix is diagonal. These simplifi- 
cations make the total covariance matrix easy to invert. 

Table 13 gives the 95% confidence level upper limits de- 
termined by MINT in the 2001 campaign. We compute 
the likelihood £ as a function of flat bandpowers 5T, 
given a uniform prior distribution. The three baseline 
lengths are analyzed separately, as the correlations be- 
tween them are in this case small enough to ignore. The 
separate likelihood functions of the mixed baselines are 
added, as their data are uncorrelated in u-v space. We 
define the upper limit so that the integral of C up to the 
limit is 95% of the integral over the whole domain. The 
mixed baselines 1540) provide the strongest limit, 
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105 /uK, on the CMB anisotropy. 

8. CONCLUSION 

The MiUimeter Interferometer has operated success- 
fully as a prototype for a high-frequency CMB instru- 
ment. MINT combines heterogeneous antennas, double- 
sideband SIS mixers, monolithic channelizers, and dig- 
ital correlators. All of these elements may be features 
in future CMB interferometers. To make a competi- 
tive detection of the CMB anisotropy, such an instru- 
ment would require wider bandwidth, longer observing 
time, and more receivers (having more distinct baselines 
would also improve the sampling in ^-space). Usable 
bandwidth in the D-ban d is ultim atel y limited by at - 
mospheric emission lines (|Matsushit'a fc Matsuo II2003|) . 
but sid eband-separating mixers— under development at 
NRAO ijKerr. Pan fc Leduc llT99^ — offer an instant fac- 
tor of two improvement, given equal system noise tem- 
peratures. Consider an improved MINT, with 8 GHz 
bandwidth (4 GHz upper and 4 GHz lower sideband) , 
four times the actual 2001 observing time, and fully 
working polarizers. Assuming 50 separate fields of view 
and the power spectrum of the 2 002 WMAP best fit 
ACDM model ( SoerEel et aP 12003(1 . the MINT window 
functions produce CMB bandpowers of 1450 and 590 
/iK^ at effective I = 900 and i = 1390. In this case, 
the improved MINT would have the raw sensitivity to 
see anisotropy at the level of 7cr and 6cr in the short 
and medium baselines. We conclude that the tech- 
niques described here constitute a promising, if challeng- 
ing, method for observing CMB anisotropics at or above 
100 GHz. 
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